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Abstract 



At finite density, charge in holographic systems can be sourced either by 
explicit matter sources in the bulk or by bulk horizons. In this paper 
we find bosonic solutions of both types, breaking a global U{1) symmetry 
in the former case and leaving it unbroken in the latter. Using a minimal 
bottom-up model we exhibit phase transitions between the two cases, under 
the infiuence of a relevant operator in the dual field theory. We also embed 
solutions and transitions of this type in M-theory, where, holding the theory 
at constant chemical potential, the cohesive phase is connected to a neutral 
phase via a 2; = 2 QCP of Schrodinger type. 



1 Introduction 



Given a holographic system at finite density with respect to some U(l) symmetry, we 
would like to answer basic questions relating to the physical nature of its groundstate 
such as 'is the symmetry broken or unbroken? ' and in the latter case 'are there Fermi 
surfaces present?' and 'what are the properties of such Fermi surfaces?'. Lately such 
questions have attracted a great deal of interest leading to general statements along 
the lines of a holographic Luttinger relation and its bosonic analogue [H El |3]. The 
Luttinger count relates the sum of the volumes of all Fermi surfaces present to the total 
charge, while its bosonic analogue relates the total charge to the 'Magnus force' felt 
by a test vortex in the dual field theory. Such studies were guided by the realisation 
that the finite charge density of a field theory - as encoded in the dual geometry 
via a finite electric fiux - can be sourced either by explicit charged matter, such as 
fermions or charged bosonic condensates, or by a charged horizon. The latter situation 
corresponds to having some of the charge density of the field theory tied up in gauge- 
variant operators, and correspondingly not contributing to a naively defined Luttinger 
count or its bosonic analogue. The charge contributing to this deficit is identified with 
so-called 'fractionalised' charge in the field theory, that is charge carried by deconfined 
'quark' degrees of freedom which are charged under the gauge symmetry. This point 
of view has recently been bolstered by the identification of the corresponding 2kF 
singularities associated with the 'hidden' fermionic degrees of freedom by Faulkner and 
Iqbal m. 

There can be interesting phase transitions between situations in which fiux is sourced 
by a horizon to ones where it originates from explicit matter sources in the bulk [H El E] • 
An early example of such a transition was observed in the M-theory superconductor [71 
[8], where it was demonstrated that the zero-temperature limit of the superconducting 
black hole was a charged domain-wall geometry without horizon due to the interplay 
of the potential and effective gauge coupling as defined by a neutral scalar field. 

The ground states of the Abelian Higgs model with a W-shaped potential were 
analysed in [9l |T0], where it was observed that the nature of the ground state can 
change quahtatively as a function of the gauge coupling. In [71 IH] it was also pointed 
out that the disappearance of the horizon as T — )■ can be related to Nernst's Law since 
the charged domain wall has zero entropy, whereas the naive extrapolation of the finite- 
temperature geometry to an extremal black hole would carry entropy. It is thus perhaps 
expected that the charged domain wall is thermodynamically preferred over the finite- 
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entropy black hole solutions in precisely that region of the phase diagram where both 
solutions coexist and thus the 'entropic singularity' of the M-theory superconductor is 
cloaked by a superconducting dome. 

In this work we explore bosonic analogues of the fractionalisation transition of [I], 
corresponding to transitions from phases with broken U(l) symmetry to phases with 
intact U(l) symmetry and all charge associated with a horizon. The former constitutes 
the analogue of the 'mesonic phas^ of [IJ and the latter corresponds to a fully fraction- 
alised phase. There may also be partially fractionalised phases, in which some of the 
charge is carried by the horizon and some by the scalar matter. Because of the presence 
of non-singlet matter, the partially fractionalised phase always breaks the U{1) symme- 
try, so that the transition between a cohesive and a partially fractionalised phase must 
occur within the superfluid state. Hence we refer to the two kinds of ordered phases in 
this paper as the 'superfluid cohesive phase' and the 'superfluid fractionalised phase'. 
We study these transitions in a variety of models, both in a 'bottom-up' and a 'top- 
down' context, with a range of different behaviours, and we describe the conditions 
under which fractionalisation transition can occur. Among other things we describe 
a supergravity-embedded superfluid cohesive to neutral phase transition, which has 
certain aspects that are reminiscent of the standard superfluid to insulator transition 
in the Bose-Hubbard model at constant chemical potential. 

It is worth commenting that the theories considered in this paper have a structural 
similarity to the Thomas-Fermi approach used in the electron-star literature [12] with 
the charged scalar field taking on the role of the fluid of bulk fermions. It should not 
come as a surprise then that some of the IR geometries encountered in those studies 
will also feature in the bosonic theories here. 

An important role is played in the present work by IR geometries with certain scaling 
or hj^erscaling properties. Such solutions have previously found applications in applied 
holography, for example in describing QCD like theories fl3[ [HI [15] , and more recently 
in describing quantum-critical matter at finite density [T6|, [T71 UHl [TH [201 [211 [22] ■ 

The paper is organised as follows. Immediately following this paragraph we briefly 
summarise our results. In section 2 we introduce the class of Einstein-dilaton-charged 
scalar models used throughout the paper and develop their thermodynamics and holo- 
graphic renormalisation. Section 3 gives a detailed discussion of certain bottom-up 
models which illustrate the general concept of bosonic fractionalisation and motivate 

^ Apart from this one instance we adopt the term 'cohesive' for any non-fractionahsed phase [TT], 
in order to avoid potential confusions with the more conventional meaning of the term 'mesonic'. 
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the more subtle M-theory case to which we turn in section 4. We finish with a discussion 
and a short appendix. 

1.1 Summary of results 

The sahent features of the analysis boil down to the IR structure of the zero-temperature 
states. We encounter three illustrated in Fig. [1} depending on the structure 

of the effective gauge coupling, which is controlled by the neutral scalar field. If the 
gauge coupling vanishes in the IR, we find a fully fractionalised phase, or a superfiuid 
fractionalised phase. If it remains finite, we find a cohesive phase. The U{1) symmetry 
is broken in the latter two phases by explicit charged matter outside the horizon. In ad- 
dition to the fully fractionalised phase, there is one further unbroken phase, where the 
neutral groundstate does not depend on the chemical potential. This 'incompressible' 
phase together with the transition into it from the cohesive phase via a, z = 2 criti- 
cal point is reminiscent of the superfiuid to insulator transition in the Bose-Hubbard 
model. 

2 General Features 

In this section we introduce the theories leading to quantum phase transitions in terms 
of a framework that is general enough to cover both the bottom-up constructions of 
section |3] and the M-theory compactifications of section |4j Our general strategy follows 
a two-pronged approach. Firstly, at zero temperature we will carry out an analysis of 
possible IR geometries and the deformations which allow them to be connected to the 
asymptotic AdS configuration universal to all solutions considered in this paper. We 
are interested both in IR geometries associated with domain wall configurations and IR 
geometries associated with (extremal) black-hole horizons. This allows us to classify 
the possible ground states of the holographic theories in this work. Secondly we will 
extrapolate the known finite-temperature solutions and their dilatonic deformations 
to very low temperatures. These can lead to a state with broken U{1) symmetry via 
a condensation instability at some critical temperature or to a state in which this 
symmetry is intact. In each case we are able to connect them to the appropriate T = 
geometry. In this way we can also study other interesting low-temperature signatures, 
such as the dynamical critical exponent z and the hyperscaling violation exponent 6, 
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class la- 



superfluid cohesive superfluid f 'ised fully f'ised 

?c ff 



fully f'ised superfluid f'ised fully f'ised 

class lb 9 

9f 9c 9f 



class // 



neutral superfluid cohesive neutral 



9c 9c 



Figure 1: Zero-temperature phase structure of the the models considered in this 
paper. Shades of red indicate broken U{1) symmetry, with dark red symbohsing 
the cohesive states and bright red the superfluid fractionalised ones. The U{1) 
symmetry is unbroken in the green fully fractionalised and grey neutral phases. 
The class lb example is even in while the class la example has no partic- 
ular discrete symmetry associated with the dilaton. Note that the M-theory 
embedding in section [4] is in class //. 



in terms of which the low-temperature entropy scales like 

S^T"^. (2.1) 



We are able to identify z and 6 from an appropriate (hyper-)scaling IR geometry, 
matching perfectly to numerical fits of our low-temperature analysis. 

2.1 Action and Field Content 

The class of theories we work with in this paper is described by an action of the form 



S = I d'^'xV^[R- \Zp{<I>)FmnF^''' - '^{V<l>f - Zs{S)\DS\' -V{<I>,\S\) 



+ 'S'cs • 

(2.2) 



We denote the d+1 bulk directions with upper-case Latin M, N, . . ., whilst field-theory 
coordinates will be denoted with lower-case Greek /i, i/, . . .. The field $ is a neutral 
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scalar field, which we will refer to as a 'dilaton', whereas is a charged scalar field 
with covariant derivative DS = {V — iqA)S. We have included Sqs, which stands for a 
possible Chern-Simons or axion-like contribution often found in consistent supergravity 
reductions. Whilst the solutions we construct do not depend on Sqs, we note that in 
certain cases the phase structure may be altered by its presence, for example through 
the spontaneous formation of stripes. 

We allow the couplings of the gauge field, as well as the charged scalar kinetic 
terms, to depend on $ and S respectively. The interesting transitions observed in the 
remainder of this paper can be traced back to the interplay of charged and neutral 
scalar couplings and the potential. Often it is convenient to put the scalar kinetic 
terms into canonical form. A little thought shows that this is only possible for the 
magnitude of the complex field. Writing in a polar decomposition with phase iqip 
and then using a field redefinition of the magnitude, we can write the kinetic term of 
the charged scalar as 

Zsi\S\)\DS\' = (Vvf + q'XivfiV^ - Af . (2.3) 

We should require that X{ri) has a small-// expansion X{ri) ~ r/ + 0{7f). The hyper- 
scaling geometries are generated by a divergent IR dilaton $ — > ±00, and so it is 
important to distinguish theories where such a divergence can happen from those were 
it cannot: 

IR 

class/: Zf{^) ■ — > 00, 

IR. 

class // : Zf{^) > (always stays bounded) . (2.4) 

Furthermore Zf{^) may or may not have even symmetry under $ — > — $ as the 
IR limit is approached, again with important implications for the phase structure at 
zero temperature. In class / we are particularly interested in Zp{^) without this 
symmetry, so that e.g. if Zf{^) diverges as $ — > +00 then it remains bounded when 
<I> diverges with the opposite sign. In this manner the dilaton can naturally drive a 
fractionalisation-type phase transition. 

Finally we require that the potential has an expansion 

£V(r/,$) = -6 -$2-2772 + ... , (2.5) 

ensuring that there exists an AdS^ vacuum solution with associated AdS length i'^, 
and that the two scalar^ $ and S are dual to A = 2 operators. For the solutions in 

^The seemingly different mass terms are caused by the different normalisations of the respective 
kinetic terms. We are using a convention that appears to be standard in the literature. 
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this paper we can take the scalar S identically real, so that without loss of generality 
(p = 0. 

Throughout this paper we use the metric ansatz 

ds' = -f(r)e-^(^Uf + + ^ {d^' + dy') ■ (2.6) 

2.2 Holographic renormalisation and asymptotic charges 

We are interested in solutions that asymptotically approach the AdS4 fixed point in 
the UV. Such solutions admit a large r expansion 

f - ^ + + + ^ + 

n \ 

(2.7) 



By evaluating the asymptotic stress tensor (which involves introducing counterterms 
to render the renormalised expressions finite) we can deduce the exact relationship of 
the integration constants introduced here and the physical charges of the dual field 
theory. This is achieved by adding the specific counterterm action 

where we have chosen the counterterms for the scalar fields such that the fixed rii, $i 
ensemble has a well defined variational principle. /C = "^^^IC^u is the trace of the 
extrinsic curvature of a constant r hypersurface with unit normal rf — f{r)^. We find 
that the total action 

S,,^^S + S,t (2.9) 
is finite, as is the renormalised stress tensor 



T — 2 



\{\S? + \^')l,.. (2.10) 
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From this we can relate the expansion coefficient Gi to the energy densitjj^ defined as 

e = ^Too = -J (Gi - r]ir]2 - |$i$2) • (2.11) 

The two scalar fields play very different roles in our system. The dilaton $ gives rise 
to an explicit deformation parameter: we switch on a relevant operator by sourcing it 
in the dual theory, i.e. by taking a non-zero value of the UV parameter $i, and then 
study the behaviour of the theory as we vary this source. In contrast, the field t] is 
only ever allowed to condense spontaneously, i.e. we set rji = 0. 



2.3 Thermodynamic Relations 

By analytically continuing to Euclidean signature 



-iSr, 



we define the grand canonical potential 

Ie = mf^.T) = Pvohuifi,T) . 
Evaluated on-shell, the action can be written as a total derivative 



5, 



OS 



i/3vol2 



d_ 

dr 



-99' 



dr 



(2.12) 



(2.13) 



(2.14) 



which vanishes at the horizon, and thus depends only on the asymptotic charges. Using 
techniques of [23l |2^ one can show that this is the negative of the pressure 



u{^i,T) = -P. 



(2.15) 



We can exploit the symmetries of the background metric to derive a Smarr-Gibbs- 
Duhem relation. This follows from the vanishing of a certain total derivativ^on shell, 
as we now show. By using the Killing equation in addition to the Ricci identity for the 
two Killing vectors dt and dy we can derive the relation 



-g {R\ - Ry) = -dr 



2P 



//?' - /' + 



2/ 



r 



(2.16) 



■^Thc quantity e"''" is the boundary speed of light (squared), and so the conversion between ADM 
mass and energy, and correspondingly energy density, involves a factor of e"^". For this and related 
reasons it is convenient to set Pa — 0, which we will assume to be the case from now on. 

*The on-shell vanishing of the total derivative explains why the authors of [3123] were able to write 
the on-shell action as two different boundary terms. 
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It is then a simple matter of using the trace-removed Einstein equations to find that 
on-shell 



1 



g {R\ - Rl) = --dr [v^Z(0)A,F*^] 



(2.17) 



so that, upon integrating from the horizon to the UV boundary, we obtain the relation 



fP' -f + 



2/ 



-gZ{<j))AF 



tr 



Evaluated on our expansions this gives the Gibbs-Duhem relation 

3 



-e = nQ + Ts 



or, written in terms of the pressure 

6 + P = fiQ + Ts . 
Here we have used that the temperature is given by 



(/3„-/3+)/2 

T = . fir) 



47r£2 



(/3a = 0) . 



r=r-|- 



Moreover, the entropy density is given in terms of the horizon radius as 



(note IGttGn = 1) • 



AG 



N 



(2.18) 



(2.19) 



(2.20) 



(2.21) 



(2.22) 



The above derivations continue to make sense for solutions of the theory (2.2) where 
lim.r_>rjj5 f'{r) = as — )■ 0, i.e. solutions of vanishing entropy at zero temperature. 
We also allow the dilaton field to diverge logarithmicalljj^ in this limit, as long as it 
does not interfere with the conditions just mentioned. 



3 Class /: A Bottom-Up Model 

When Z is not symmetric in $ there is the possibility of two qualitatively different 
behaviours assosciated with a dilaton that diverges logarithmically to either plus or 
minus infinity. In the following models, we will see that in the first case the effective 
gauge coupling ~ Z^^ will go to zero in the IR and in the second it will diverge. There 
is a continuous transition between a superfluid fractionalised phase and a cohesive 
phase, via a quantum critical point with a finite dilaton. Both phases will be of 

^Specifically the dilaton blows up logarithmically as a function of radial proper distance. 
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hyperscaling form. This phase transition will be driven by the asymptotic value of the 
dilaton. We will also see an additional phase transition between a partially and a fully 
fractionalised phase, assosciated with the edge of the superconducting dome. 
We will consider the class of models 

= Zo'e-^*/^ , fV{^, ri) = -V^^ cosh (fe^/v^) - 2ri^ + gy , (3.1) 

where a,b > 0. The rj'^ term allows for interpolating soliton solutions between the AdS^ 
maximum at r/ = and the AdS^ minimum at rj"^ = g'"^. We will study IR geometries 
for general values of a and b (see also [U E]). All our numerical results will be for the 
particular case a = b = g^ = 1, Vq = 6 and Zq = 1. 

3.1 Zero Temperature 

We now study the phases of our system at zero temperature. We construct the IR 
geometries by making a scaling ansatz for the fields, allowing for a possible logarithmic 
divergence of the dilaton. With this divergence the potential is schematically V ~ 

_l_ ^-s^ ^Qj. gome 6, and so these scaling solutions are not exact: they are the leading 
order parts of series solutions, with subleading corrections in increasing powers of r. 
For each of our three cases we also identify the possible irrelevant deformations. 

The divergent dilaton endows our leading order solutions with hyperscaling symme- 
try [201 EI], which is a generalisation of the familiar Lifshitz scaling: 



X — 

t 

r — i- 

under which ds — )■ X^^'^ds. Here z is the 
the hyperscaling violation parameter. 

3.1.1 Fractionalised IR Solutions 



Ax , 

A(^-2)/V, (3.2) 
usual dynamical critical exponent, and 6 is 



In this case some or all of the flux is sourced by the horizon. Starting with the fully 
fractionalised case, we switch off the condensate. We find the following scaling solutions 
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at leading order: 



with 



fo 



2(a-h) 

f{r) = for , 
At{r) = Aor , 

$(r) = ^logr, (3.3) 

a + 



4(6 + a{a + 6))(12 + (3a - 6)(a + 6)) 



^ 4(6-&(a + &)) 

° (12 + (3a-6)(a + 6))Z2- ^ ' 

These solutions have dynamical critical exponent z and hyperscaling parameter 6 given 

^ ^ 12 + (a-3|a + t) _ ^ «^ 

— 0^ b — a 

Note that though these parameters are infinite in our special case a = b = 1, their 
ratio is still finite (such a situation is discussed further in [H]). This ensures that for 
solutions close to this T = geometry, the various thermodynamic parameters have 
finite scaling exponents with temperature. 

There are two irrelevant deformations about this solution. The first simply corre- 
sponds to a shift in the IR vacuum expectation value of the charged scalar: 

Sv = Vo- (3.6) 

Integrating this deformation to the UV gives geometries with explicit bulk charge in 
addition to horizon charge. We thus find that the fully fractionalised solutions connect 
continuously to the superfiuid fractionalised phase branch of solutions. The second 
deformation is assosciated with the dilaton, and has the form 

5$ = a|,^r^*, 5/ = aJ«r^/, 5/3 = aJ^r^*, SAt = a^/r''^, (3.7) 

with exponents 

2(a - b) 



Vf = + 



a + b 



_ 12 + {3a -b){a + b) 
i^A - + , (3.8) 
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where the coefficients aj and the exponent z/$ are consistently determined by the IR 
hmit of the hnearised equations of motion. For example, in the simplified case a = b 
the theory admits an irrelevant deformation with exponent 

3 + 6^ 

and coefficients 



(^-3 + V81 - 2462 j ^ (3 9) 



{ttf , , , a$} = a^^ |-/o , 2 , ^^'^^^2^ ' "Tf} ' ^^'^^^ 

with an overall free magnitude a^^. This is to be compared with the analogous mode 
in p. 

3.1.2 Critical IR Solution 

The critical IR solution has non-divergent dilaton (in fact $ = 0) and is just the AdSi 
minimum that arises when the condensate sits at the bottom of its potential, rj = ±(7^ ^, 
with vanishing gauge field At{r) = 0. Consequently its AdS length scale is shifted away 
from the usual value by a (7,?- dependent factor: 



/(r)e-^(") = 
fir) 



^2 



Vir) = ±9-\ (3.11) 

with 

There are two irrelevant deformations about this solution: one involving the flux and 
one the charged condensate. They have the form: 

5 At = afr""^ , 5r] = a^f^r''^ , (3.13) 

with exponents 



^i-i + .i' '''''' 



U V ^0(1+9?^.'), 
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with 



3.1.3 Cohesive IR Solution 

We seek cohesive solutions by demanding that the horizon flux limr_s>o a/— 
vanishes. At leading order we flnd the following scaling solutions and their deforma- 
tions: 

2(3-6^) 

/(r) = /or 3 , 

Atir) = 0, 
2b 

$(r) = — logr, (3.15) 
v3 

31/2 

These solutions have dynamical critical exponent z and hyperscaling parameter 9 given 

by 

^ = 1, () = -J^,- (3.17) 

Just as in the fractionalised case, the vacuum expectation value of the charged scalar 
is not flxed, St] = tjq. There is then a second irrelevant deformation which introduces 
the required flux. For a = b this is 

Mi = r'^+^^('^o^ (3.18) 

This exponent depends on the IR value of the charged scalar: 



This is merely the leading order contribution to the series solution from the deforma- 
tion; there are additional subleading corrections in powers of 21/ a- We thus have the 
non-trivial consistency condition that this power be positive, which translates into the 
constraint on the condensate: 



2 (3 + 62)2/72 

Vl > ^ J/ ° • (3.20) 
A similar structure was seen in the zero-temperature superconducting solutions of |25j . 
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3.2 Finite Temperature 

In order to complement our zero-temperature picture, and to complete our phase dia- 
gram, we would like to study how these solutions extend to finite temperature. Holo- 
graphically, a non-zero temperature is indicated by the presence of a (non-degenerate) 
bulk horizon. We thus make the ansatz 

/ = f+{r - r+) + ■ ■ ■ , 

/3 = /3+ + ■ ■ ■ , 

At = A,+ (r-r+) + -- - , 

$ = $+ + ■■■, 

V = V+ + ---- (3.21) 



Temperature masks the singularities of the solutions of section 3.1, and the dilaton 
is finite here. Despite this we will see clearly the imprint of the zero temperature 
IR scaling behaviour at finite temperature. In particular, our results will exhibit the 
striking infiuence that the quantum critical point has on the finite temperature physics. 

Such solutions exhibit a superconducting instability I2SIEZ1, whose critical temper- 
ature depends on the value of the dilaton deformation This critical temperature 
can be dialled to zero, forming the edge of a superconducting dome. We have checked 
that the free energy of the broken phase solution is lower than that of the unbroken 
phase whenever it exists. 

Since we are holding the theory at finite chemical potential, it is natural to consider 
the phase diagram as a function of $i/yU. We can illustrate the cohesive to partially 
fractionalised critical behaviour at finite T by approaching T = from the top of the 
superconducting dome (represented by a black dashed line in each case) at constant 
values of $i/yU. A selection of such slices is shown in Fig. [2] We see that the critical 
solution approaches the unique constant dilaton solution in the IR, while solutions on 
either side of it exhibit an IR divergence as T = is approached, consistent with our 
earlier analysis. 

A comprehensive visualisation of the critical behaviour is given in Fig. |3] where we 
show a fit to the scaling a of the entropy density of the form 

4 ~ ci (-V , (3.22) 



for the region below the superconducting dome. From equation 2.1 we expect that near 
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Figure 2: Cohesive to partially fractionalised critical behaviour in the broken 
phase. The dashed black line indicates the edge of the superconducting dome 
and the coloured lines are some slices at fixed $1/// near the critical solution 
(black). The black dot shows the critical T = solution with good agreement 
with the critical T > solution. 



T = 0: 

a ~ . (3.23) 

The emergence of the quantum critical wedge associated with the z = 1 quantum 
critical point of the cohesive to partially fractionalised ground state is clearly visible. 
On either side of the wedge, the values of the coefficient a are in agreement with 
the predicted values of a = 2/3 in the cohesive phase and a = 2 in the (partially) 
fractionalised phases. 

3.3 The Fractionalisation Transition 

Our analysis thus far can be taken as strong evidence for the existence of a T = phase 
transition between a superfluid cohesive phase and a superfluid fractionalised phase at 
a critical UV deformation ^i/fi{gm) — —0.131. In this section we complete the physical 
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0.75 



Figure 3: Scaling behaviour for the quantum-critical region associated with the 
z = 1,6 = QCP. Colour indicates the value of a, which is proportional to the 
logarithmic derivative of entropy with respect to temperature at fixed /i. The 
T = QCP is located at the black dot. Hatching indicates low temperature 
regions not covered by our data. 



picture, deforming the IR geometries of section 3A by irrelevant deformations in order 
to construct full solutions which asymptote to AdS4 in the UV. 

Our numerical results are conclusive: we find one-parameter families of full ge- 
ometries of each of the three types — superfluid cohesive phase, superfluid fraction- 
alised phase and fully fractionalised. The superfluid cohesive phase exists for < 
^i/f^idm)- The superfluid fractionalised phase exists for < < ^i/fi{gf), 

meeting the fully fractionalised branch smoothly at 0.621. 

A holographic measure of fractionalisation is the ratio of flux emanating from the 
deep IR of the solution {e.g. a black-hole horizon) and the total charge Q: 



A ^ 1_ 
Q ~ Q 



r+ 



gZ{^)F''' . (3.24) 



As a consequence of Gauss's Law, this must interpolate between zero for the superfluid 
cohesive phase and unity in the fully fractionalised phase [H [3l [21] . In Fig. Ill we plot 



15 



this measure for each of the solution branches found, which confirms their identities as 
either superfluid cohesive, superfluid fractionahsed or fully fractionalised. 




Figure 4: Dependence of the T = domain walls on the parameter As 
in Fig[l| shades of red indicate broken U{1) symmetry, with dark red indicating 
a superfluid cohesive and bright red a superfluid fractionalised phase. Fully 
fractionalised geometries are shown in green. Supplementary low temperature 
data (T ~ 10~^/i) in the vicinity of the fractionalisation transition gc is indicated 
by black dashes. 

For the above phase transitions to occur, it is necessary that the fully fractionalised 
branch be thermodynamically subdominant wherever it coexists with the superfluid 
cohesive or superfluid fractionalised branches. To this end, we plot the free energy 
density u of each family in Fig. |4j We see that the free energy is indeed higher in the 
fully fractionalised phase. It proved numerically unprofitable to calculate the T = 
free energies of the cohesive and superfluid fractionalised branches in a small region 
close to the critical point, so in that region we have supplemented our results with low 
temperature (T ~ 10~^/i) data. Furthermore we note the precise agreement between 
the T = analysis and the low temperature results, which is thanks to the plateauing 
of the thermodynamic quantities near T = that we (numerically) observed in Fig. [2] 

3.4 Class lb: A Bottom-Up Model 

The results we have presented so far map out the detailed phenomenology for a theory 
in class la, under the classification of Fig. [ij for which the coupling term Zp remains 
finite for an interval of UV scalar deformations $i. Consequently, we argue, the theory 
in this interval was able to support superfluid cohesive solutions which provided the 
thermodynamically preferred solutions. 

Here we make a few brief comments on a model governed by the same potential 
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(3.1), but with a $-even gauge coupling term Zp{^) = Zq cosh with a > 0. The 
important feature of this model is that Zp will diverge in the IR given an IR divergent 
$, irrespective of its sign. Thus, the interval of cohesive solutions exhibited by our 
earlier model will here be reduced to a single point, $i = 0. At this point the bulk 
solution has $(r) = everywhere and is simply an AdS4 to AdS4 charged domain wall. 

Thus our model exhibits a transition from superfluid cohesive phase at the point 
$1 = 0, to superfluid fractionalised phase < |$i| < ^i{gf), and ultimately to a fully 
fractionalised phase when |$i| > ^i(gf). This is illustrated in Fig[l| We now turn to a 
model which arises from a consistent truncation of IID supergravity, which shares even 
gauge coupling Z with class lb, but which differs from all models in class I by having 
finite Z everywhere. As we shall see this leads to markedly different phenomenology. 



4 M-Theory 



We work with the consistent truncation of [HI EH] , whose equations of motion can be 
obtained from the action 



S 



M 



WttG 



2(l-/l2)2 



R- 



1 + 3/^2 



77 T^MN 



2(l-||xlT 



\Dx? 



6 {-l + h^ + 


X 




£2(1-1 


X 


2)2(1. 


- /l2)3/2 



+ 



IGttG 



2h{3 + h'^ 
l + 3/i2 



-F AF. 



(4.1) 



By defining x = ^e^'^'^ and performing the transformation 



h = tanh 



p = —;= tanh ^ ^ 



this action falls into the class II of models of ( 2.2[ 2.3) with the specific choices 

4 

V2 



(4.2) 



Zpi^) 



-^h' (A) (l + 3tanh^ (^) 



X{r]) = —= sinh (^V^f] 



(4.3) 



1, This model also has a Chern-Simons term. We now turn to 



and charge qi 

a description of the detailed phase structure of the CFT duals of this system, held at 
finite density, with particular regard to the zero-temperature ground states. Much of 
this discussion has been reported in the literature [H [8] , and more details can be found 
there. 
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4.1 Phases 



As before, a superfluid branch of solutions emerges as an instability of the charged 
Reissner- Nordstrom family of solutions which exists in this model when h = 0. Once 
again, under the dilaton deformation in the constant fi phase diagram we see a dome of 
superconducting solutions, whose T = limit, as shown in [TJ E], is a charged domain 
wall interpolating between two AdS regions of different AdS lengths. Whenever the 
U{1) symmetry remains unbroken the zero-temperature limit of the neutral and charged 
black hole solutions is incompressible, in the sense that it does not depend on the value 
of fi. In fact it approaches a zero-density state with vanishing (direct) conductivity. 
Furthermore, the T — )■ limit of the neutral and charged solutions meet at the unique, 
unbroken T = solution of the system, which lifts to an eleven-dimensional Schrodinger 
solution of M-theory [29j. Thus the entire regioi]|^ at T = outside the dome is 
degenerate on a constant-/! phase diagram. In view of this 'incompressibility' of the 
groundstate. Fig. |5] shows the domains where these solutions exist at fixed energy 
rather than fixed chemical potential as in [S]. 

As we tune the deformation by Oh, both at T = and at T 7^ 0, the expectation 
from our previous results in this paper is that we eventually reach a transition from 
the broken U{1) solution to a (partially) fractionalised phase. 



We now turn to an exploration of the neutral solutions of (4.1) outside the dome in 



which the neutral scalar h attains the singular value h = 1 in the IR. In the parametri- 



sation of Eq. (4.2) we see that this corresponds to a (logarithmically) divergent dilaton 



singularity, which should by now be no surprise. 
4.2 Neutral top-down solutions 



We rewrite (4.1) for the case of interest assuming the charged scalar x is trivial, since 



we do not want solutions in which the U (1) symmetry is broken. Hence, passing to the 



variables (4.2), we find 

/ d^xV^glR - """^ (f/ - hv^r + I cosh U)] . (4.4) 

J L l + 3tanh2(^^j 2 P Wsyj 



S 



Recall that at finite temperature we have a one-parameter family of neutral 'dilatonic' 



black holes, with IR behaviour (3.21), where ^ = 0. At zero temperature the nature of 



^The possible emergence of non-isotropic phases [30l|3T] is beyond the scope of this work, but could 
modify the picture in some regions. 
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& charged domain walls 

Figure 5: The existence of bulk solutions at fixed mass. The red line denotes 
neutral solutions terminating in the Schr zero-termperature fixed point. Note 
that the neutral dome meets the superconducting dome at precisely this fixed 
point, showing that there is no fractionalised phase in this model. 



the solution changes, as we now describe. Guided by our findings above we suppos^ 
that the isometries of the metric are enhanced from the M x 5*0(2) symmetry above 
to the full S0{1, 2) Lorentz symmetry. By choosing a convenient radial gauge we take 

ds"^ = + p^sech (<!?{p)/V^) v^.udx^dx" . (4.5) 
F{p) \ J 

This expression is convenient for the current purposes, but we note that we can convert 



it into the form (2.6) via 

dp 



p'sech^^/Ts) , e-^f = r\ (4.6) 



We can use this ansatz to construct the singular IR solutions admitted by the M-theory 
truncation. 



''In support of this assumption we have strong numerical evidence for this behaviour by studying 
the zero-temperature approach of the none-extremal black hole solutions. 
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4.2.1 Singular IR solutions 



We can construct the desired IR solution by a standard dimensional reduction trick. 
Suppose then that $ diverge^ in the IR as $ — > oo. In this limit the action approaches 



smg 



L 2 . 



which can be obtained in a reduction via a gravi-photon ansatz 



(4.7) 



(4i 



where A G T*M and hats denote five-dimensional quantities. The starting five- 
dimensional action is simply Einstein-Hilbert with a cosmological constant: 

3 



d^x 



R + 



(4.9) 



Under this reduction, pure five- dimensional AdS with metric 

'dp' 



dsl = il 



+ p' {-dt' + dx') + p'dz' 



(4.10) 



reduces to a logarithmic dilaton solution in four dimensions, that is: 

^dp' 



ds' 



P 



+ p'^-dr + dx.' 



with $ = -V3\og{p) 



(4.11) 



This solution, as we shall now see, plays the role of the singular IR we are after. Note 
that in this case the gravi-photon field A is trivial, so that the IR solution is obtained 
via a simple circle reduction from five dimensions. 

Note that we also have a five-dimensional Schrodinger metric [321 [331 EIIj obtained 
as a simple deformation of the above, by adding a term proportional to the light-cone 
coordinate {dx~^)': 



ds' 



-P'p\dx+)' + ^+p' idxl + 2dx+dx-) 
pi 



(4.12) 



where /3 is a constant. Once reduced to four dimensions and seen as an expansion 
in small p, the Schrodinger and AdS metrics in fact agree to leading order, that is 
they give the same IR behaviour for the desired logarithmic dilaton solution in four 



dimensions. Of course, solving the full equations coming from (4.1 ) order by order will 



^The case $ — >■ —oo is equivalent by the $ — > — $ symmetry of the lagrangian coming from IID. 
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determine the exact solution and thus distinguish between the two cases. This was 
pointed out earher in [29] . 

Here we demonstrate that this solution is in fact the universal T — )■ limit of the 



model (4.1) outside the dome. The full M-theory solution, lifted to eleven dimensions, 
is a Schrodinger solution Schrs x KE^ with dynamical critical exponent -2 = 2, where 
the Schrs is fibred non-trivially over the Kahler- Einstein space KEq. In the dual field 
theory this solution corresponds to a mass-deformation driving the system to the non- 
relativistic z = 2 YR fixed point [29J. The full structure of this is not evident from 
the four-dimensional perspective, but we see the z = 2 scaling for example in the 
temperature scaling of entropy density at criticality. 

Returning thus to constructing the neutral IR solutions, we substitute the metric 



ansatz into the equations following from (4.1), resulting in a single, fully decoupled 



non-linear ODE for the scalar field $. Its precise form is given in appendix |Aj Once $ 
is determined, the Einstein equations determine the metric function F[p) algebraically: 



. ^ 4rVcosh($/v/3) 

4-(p$VV3)2sech2($/\/3) -4p<l>'/\/3tanh(<l>/v^) ' ^' ' 



where the prime denotes a derivative with respect to the radial direction p. The $ 



equation admits an IR expansion, which reproduces (4.11) to leading order. It takes 
the form 

$(p) = -v/31ogp+^p2 + ---, (4.14) 

and contains no free parameters. Note that the freedom associated with the p rescaling 
of the metric has no physical effect on the solution. This unique IR expansion can be 
integrated to the UV and connects to the desired skew-whiffed AdS^ solution. We 
have done this and read off the UV data. In appropriate unit one finds h2/h\ = 



-0.288... l/(2v^) for the expectation value (Oh) and e/hl = 0.779... ~ ^ 

for the energy. The analytic values are those of the exact solution found in [S]. One 
can check that the chemical potential p{r-oo) — p{tyb) vanishes. We thus conclude that 
the unique neutral T = solution of the theory ( |4.1 ) with logarithmically diverging 



dilaton is the analytic solution given in Eq (8.2) of |H]. 
4.3 Approaching the critical point 

It is enlightening to study the power-law associated with the approach to zero tem- 
perature. At the critical value we expect scaling behaviour compatible with the z = 2 
^We quote numerical values in the same convention as [Sj for ease of comparison. 
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fixed point, in other words T-linear scaling. Consulting Fig. |6]we see that this scaling 




Figure 6: Approaching the M-theory Schrodinger critical point from finite 
temperature. The black line for /ii//U = /iic/^ = 0.354... shows precisely the 
critical T-linear scaling. It delineates the red supercritical line {hi/fi = 0.382) 
from the blue subcritical line {hi/ ^ = 0.324). Note that the complex scalar field 
is set to zero in all cases here, and the subcritical region is masked in the full 
theory by a superconducting dome. 

behaviour is indeed true precisely, if the temperature is lowered holding at the 

critical value. The critical curve in black clearly approaches z = 2 m the appropriate 
limit. The red curve corresponds to the approach of the (unstable) AdS2 infrared, and 
so the diverging 2; — )■ oo is expected. Note that this region is cloaked by the super- 
conducting dome in the full phase diagram. The blue curve corresponds to a value 
of hi/^ outside the dome. The value of a remains bounded as T approaches zero in 
this case. While the zero-temperature solution itself, for hi/n> hid n is independent 
of the value of the chemical potential the finite-temperature solutions are not, and 
so the approach does depend on ratios like hi/ or T/yU. This is clearly illustrated 
in Fig. |6| where the critical scaling with a = 2 delineates a region of subcritical (as 
exemplified by the blue curve) scaling with a <2 from a region of supercritical scaling 
with a — )■ 00, as exemplified by the red curve. The latter behaviour corresponds to 
approaching the zero-temperature AdS2 geometry, which in the full theory is masked 
by the superconducting dome. 
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5 Discussion 



The results of this paper are twofold. On the one hand we have constructed bosonic 
analogues of so-called 'fractionalisation transitions' by relying on suitably constructed 
minimal bottom-up models. On the other hand we reexamined the phase diagram of 
the M-theory superconductor in this light, showing how the analogous transition in 
this system is in fact a transition from a superfluid cohesive phase to a neutral phase 



with vanishing (direct) conductivity We emphasise again, that these transitions are 
taking place at constant chemical potential. 

The IR behaviour of the dilaton coupling to the gauge kinetic term at zero tem- 
perature governs the quantum phase transitions observed in this paper. All examples 
encountered in this work support the conjecture that one cannot have a cohesive phase 
if the dilaton coupling to the gauge-kinetic term diverges {Z — ?■ oo). While we do not 
have a rigorous argument that would prove this, we can see heuristically that it makes 
sense: diverging Z means that the effective gauge coupling vanishes as is obvious in a 
normalisation where the gauge-kinetic terms reads ~ ^F^. That means that matter in 
the strict IR, where Z diverges, cannot source any flux. Thus any flux there can only 
emanate from a horizon and consequently we would have a fractionalised contribution 
to the overall flux. 

In asymptotically AdS spacetimes, one often finds that the zero temperature limit 
of a black hole solution has an emergent AdS2 geometry, indicated by a divergent 
dynamical critical exponent. The dual field theory then has a non-zero entropy at 
zero temperature. Often these finite entropy geometries are unstable when embedded 
in top-down models (see e.g. [71 El ED]), so that the third law of thermodynamics is 
upheld by the classical geometry. 



As is clear from equation (2.1), and strikingly visualised in Fig. |3| the T = 
infared hyperscaling parameter governs the scaling behaviour of various thermodynamic 
quantities at low temperature. In particular, a divergent critical exponent can be 
compensated for by an equally divergent hyperscaling parameter [TT], resulting in a 
vanishing entropy at T = 0. In this manner one can avoid instabilities associated 
with finite entropy. Our bottom-up model of section [3] uses precisely this mechanism. 
The M-theory model, however, is of class II, and so there is no fractionalised phase 
associated with a vanishing gauge coupling and its corresponding hyperscaling violation 
solutions; we indeed see that the fractionalised (Reisner-Nordstrom) phase is unstable 



^^Interestingly, however, there is a non-vanishing transverse conductivity Heaxy 0. 
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to scalar condensation, and so masked by a superconducting dome. It would be very 
interesting to investigate precisely what conditions have to be met for any entropic 
singularity {i.e. any AdS2 IR geometry) to be unstable [33] to new phases cloaking it 
at low temperatures. 

Recently, following earlier work on Lifshitz geometries [36], it has been suggested 
[371 138] that hyperscaling IR geometries are destabihsed due to the running dilaton. 
If this is the case for the geometries considered here then our results will still apply 
for an intermediate range of energy scales, above the scale where quantum corrections 
smooth out the IR geometry. It is interesting in this context that at least some of these 
seemingly singular geometries can be lifted to perfectly regular higher-dimensional 
geometries, potentially eliminating the need to consider the corrections of 155] . 
Note, however, that in the case of Lifshitz geometries, the higher-dimensional lifts were 
found to be as problematic as their lower-dimensional descendants |39] , as manifested in 
the propagation of test strings. The status of the possible resolution of such singularities 
(e.g. by matter sources |1D]) in string theory deserves further attention. 
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A The ODE in section 4.2.1 



Here we give the details of the ODE determining $(p) in section 4.2.1, which follows 



from the lagrangian (|4.7|) upon setting At = jj with fx constant. 



12- p'^'ipY 4 + 3sech 



V3 



1 + 3sech^ (^)] p$'(p) + Isech^ (^) p'^''{p) . (A.l) 



For any solution the chemical potential in the boundary theory is arbitrary, as A is 
pure gauge in the bulk. If A is to be manifestly well-defined in Kruskal coordinates 
one should choose the gauge At = 0. 
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